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We use the gradient flow for the renormalization of the Polyakov loop in various representations. 
Using 2+1 flavor QCD with highly improved staggered quarks and lattices with temporal extents of 
Nt = 6, 8, 10 and 12 we calculate the renormalized Polyakov loop in many representations including 
fundamental, sextet, adjoint, decuplet, 15-plet, 24-plet and 27-plet. This approach allows for the 
calculations of the renormalized Polyakov loops over a large temperature range from T = 116 MeV 
up to T = 815 MeV, with small errors not only for the Polyakov loop in fundamental representation, 
but also for the Polyakov loops in higher representations. We compare our results with standard 
renormalization schemes and discuss the Casimir scaling of the Polyakov loops. 


I. INTRODUCTORY REMARKS 

The Polyakov loop in fundamental representation is an 
order parameter of deconfinement in SU(N) gauge theo¬ 
ries. For the SU(3) gauge group it is dehned as 

i 3 (x) = ^TrT’exp T)t°‘dr^ , (1) 

where x is the spatial coordinate, V is the path ordering 
operator, and the Euclidean time r is integrated up to 
the inverse temperature. The nonzero expectation value 
of T 3 (x) above the transition temperature signals de¬ 
confinement and screening of static quarks. After proper 
renormalization the logarithm of the Polyakov loop gives 
the free energy of a static quark in temperature units 
mm- In the confining phase below Tc the correspond¬ 
ing free energy is infinite. Above that temperature it 
becomes finite due to color screening. 

In QCD the Polyakov loop is not an order parameter, 
its expectation value is nonzero at any temperature as 
static quarks can be screened by dynamical quarks al¬ 
ready in the vacuum, i.e. the free energy of the static 
quark is always finite. Nonetheless, the temperature de¬ 
pendence of the Polyakov loop reflects the change of the 
screening properties in the medium and thus is linked to 
deconfinement. 

So far, we discussed the Polyakov loop in the funda¬ 
mental representation. One can define the Polyakov loop 
Ln (x) in any representation n by replacing the generators 

of the fundamental representation by the generators of 
the corresponding representation , as well as the corre¬ 
sponding normalization of the trace in Eq. Q, and con¬ 
sider the free energy of the color charge in representation 
n. The color charges in higher representations may be 
screened at any temperature already in pure gauge the¬ 
ory. However, also in this case the temperature depen¬ 
dence of (L„(x)) , or equivalently of the corresponding 
free energy F^, is sensitive to the screening properties of 
the medium and thus to deconfinement. 

As stated above, the expectation value of the Polyakov 
loop Pn{T) = {Ln{^)) requires renormalization in order 
to be interpreted as the free energy of static charges. The 


renormalization of the Polyakov loop is multiplicative [3] 

P„(T) = P,^“(T) = , 

( 2 ) 

where is the bare gauge coupling corresponding to a 
given lattice spacing[2]. In the fundamental representa¬ 
tion the renormalization is usually achieved by requiring 
that the free energy of a static quark antiquark pair is 
equal to the corresponding zero temperature potential 
at very short distances and assuming a certain normal¬ 
ization of the zero temperature potential. The constant 
63 ( 50 ) in Eq. (§ corresponds to the additive shift of the 
zero temperature potential ensuring that it has the pre¬ 
scribed value in physical units. Thus the calculation of 
the renormalized Polyakov loop requires the calculation 
of the zero temperature potential for each value of the 
bare gauge coupling 50 used in finite temperature calcu¬ 
lations. For higher representations one can proceed in 
a similar manner to obtain the renormalization constant 
Zn but usually, as we discuss later, the assumption of 
Casimir scaling is used to estimate them. 

The renormalized Polyakov loop in the fundamental 
representation has been calculated in SU(N) gauge theo¬ 
ries m HHS] as well as in QCD with two and three quark 
flavors with relatively large quark masses [3 [g. Re¬ 
sults for the renormalized Polyakov loop also exist for the 
physically relevant case of 2+1 flavor QCD with physical 
or nearly physical quark masses PIT5] . The Polyakov 
loop in higher representations has also been studied in 
pure gauge theory mm and in two-flavor QCD with rel¬ 
atively large quark masses 0. 

In this paper we calculate the renormalized Polyakov 
loop in 2+1 flavor QCD with physical quark masses in 
various representations. Calculations will be performed 
at several lattice spacings in order to control the dis¬ 
cretization effects. A new method for calculating the 
renormalized Polyakov loops based on the gradient flow 
[16] is introduced. The rest of the paper is organized as 
follows: In the next section we discuss the lattice setup. 
In Sec. Ill we discuss the renormalization of the Polyakov 
loop in the fundamental representation using the gradi¬ 
ent flow. The Polyakov loop in higher representations 
and Casimir scaling is studied in Sec. IV. Finally, Sec. 
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V contains our conclusions and some technical aspects of 
the calculations are presented in the Appendices. 


II. LATTICE SETUP 

We perform lattice calculations in 2+1 flavor QCD us¬ 
ing highly improved staggered quarks (HISQ) [T7] with 
lattice sizes of x At-=24^ x 6, 32^ x 8, 40^ x 10 
and 48^ x 12 and gauge configurations generated by the 
HotQCD collaboration using a physical strange quark 
mass rris and degenerate up and down quarks with masses 
mu = md = mi = TOs/20 [HIITS] . This setup corresponds 
to a pion mass of 161 MeV in the continuum limit [14]. 
The temperature T = l/(a7VT-) is varied using the lattice 
spacing a and the lattice spacing itself has been deter¬ 
mined using the ri scale defined in terms of the static 
quark potential 


dr 


= 1 


(3) 


We use the parametrization of a/ri from Ref. [T 8 | and 
ri = 0.3106 fm to convert to physical units. We will 
cover a temperature range of T = 116 MeV up to 
T — 815 MeV. In the low temperature region and in 
the transition region we also perform calculation using 
mi = ms/40 and the HotQCD gauge configurations from 
Ref. [13]. 

On the lattice the local Polyakov loop T(x) is given by 
the traced product of all temporal links 1/4 (x,r) at the 
spatial point x 


1 

TrJ]C^ 4 (x,T), (4) 

T = 1 

with 1/4 (x, t) G SU(3) and W the temporal extent of the 
lattice. Here L 3 “®(x) denotes the unrenormalized (bare) 
Polyakov loop in fundamental (“3”) representation. One 
usually considers the spatial average when calculating the 
expectation value of the Polyakov loop, 



We consider Polyakov loops L„ in higher representations, 
n = 6 , 8 , 10, 15, 15', 24 and 27. These can be constructed 
from the Polyakov loop in the fundamental representa¬ 
tion using group theory as follows [S]: 


Ai5' 

d^24 

L 27 


T^(^3 ■ ^10 ~ ^15) 5 

15 

( 10 ) 

— (^3 • ho - Iq) , 

( 11 ) 

^( ^ 6 ^ - ^8 - 1 ) , 

( 12 ) 


where I 3 = 3 + 3 “® and ^ is its complex conjugate. From 
the above equation it is clear that the Polyakov loops 
in all representations are normalized by the dimension 
of the representation and thus will approach one at very 
high temperatures. 


III. GRADIENT FLOW RENORMALIZATION 

As discussed above, the Polyakov loop needs multi¬ 
plicative renormalization. We use the gradient flow to 
renormalize the expectation values of Polyakov loops. 
The gradient flow is defined by the differential equation 
(equation of motion) [TB] 

Vt{x,n) = -gl{d^^^S[Vt])Vt{x,y) , (13) 

where go is the bare gauge coupling and S'[Vt] is the Yang- 
Mills action. The field variables Vt[x = (x, t),^) are 
defined on the four dimensional lattice and satisfy the 
initial condition 


Vt{x,g)\t=o = U^{x) , 


(14) 


with U^{x = (x, r)) being the usual SU(3) link variable 
and < is a new index for the evolution in flow time and 
has dimension [a^]. So far we have not specf.ed the dis¬ 
cretization scheme for the Yang-Mills action. We could 
use the simple Wilson gauge action m or the tree-level 
improved Symanzik action [ 20 j for One usually 

refers to these schemes as Wilson- or Symanzik flow. In 
our study we use the Symanzik flow, i.e. S')!/)] is the 
tree-level improved Symanzik gauge action. We also per¬ 
formed some calculations using Wilson flow, which are 
discussed in Appendix [^ The differential equation is 
solved using a Runge-Kutta like scheme up to the desired 
value of t. For the flow evolution of the gauge configura¬ 
tions and for the calculation of the observable we use the 
MILC code EH. 


Since Eq. (131 has the form of a diffusion equation. 


the gradient flow smears the original field Ufj,{x) at the 
length scale 


/ = . 


(15) 


Lio = • h — h) , 

Lib = Y^(^3 ■ ^6 — h) 1 


( 6 ) 

(7) 

( 8 ) 
(9) 


For this reason operators evaluated at nonzero flow time, 
i.e. operators that are constructed from Vt(a;,/x)|t>o in¬ 
stead of U{x,g) do not require renormalization El], the 
short distance singularities are removed by the gradient 
flow. Furthermore, it can be shown that renormalized 
operators at t = 0 are equal to operators at nonzero flow 
time up to multiplicative constant EU E3] if the flow 
time is sufficiently small, a <C / <C Aqcd- Therefore, 
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Figure 1. The unrenormalized fundamental Polyakov loop (left) and the renormalized Polyakov loop corresponding to flow 
time f = fo (right) as a function of the temperature T for all lattice ensembles. 



Figure 2. The free energy F 3 obtained from the gradient 
flow compared to the continuum results for F 3 in conventional 
renormalization scheme from Ref. m- 


the renormalization of the Polyakov loop can be achieved 
by replacing the original link variables in Eq. with 
Vt{x, fj.)\t>o- The choice of flow time / = -s/St corre¬ 
sponds to a particular renormalization scheme as long 
as 


a< /< l/T . (16) 

To demonstrate the above point in Fig. we show the 
bare and the renormalized Polyakov loop constructed 
from Vt{x^)\t>o corresponding to the choice / = /o = 
0.2129 fmR In what follows we will give the flow time in 


^ The value of /o corresponds to the lattice spacing of the 32® X 8 
lattices at the lowest temperature used in our analysis, and thus 
it provides a natural unit for the flow time /. 


units of /q. One can see that the strong dependence 
of the bare loop is gone in the renormalized Polyakov 
loop as expected. 

For the calculation of the renormalized Polyakov loop 
in an extended temperature region we need to change / 
such that the constraint given by Eq. (161 is always sat¬ 
isfied. To do so we proceed as follows: We define regions 
where the flow time / is constant in physical units, which 
means that changing the temperature T via the lattice 
spacing means changing the flow time t in the actual cal¬ 
culation such that / = const in fm. Different choices of / 
correspond to different renormalization schemes. For this 
reason the free energy should be independent of the flow 
time up to a constant shift, i.e. F„(/) — F’„(/') is approx¬ 
imately T-independent. We are limited in the range of 
Nr and therefore, as we want to cover a very broad tem¬ 
perature range, we have to define different flow regions 
to fulfill this condition: 


( 3/0 for T < 200 MeV , 

I 2/0 for 200 MeV < T < 300 MeV , , , 

] O. 5 O /0 for 300 MeV < T < 600 MeV , 

I O. 25/0 for T > 600 MeV , 


where /o = 0.2129 fm. The different regions can then be 
matched by a constant shift of the free energy and we 
do this by determining the shift via an overlapping tem¬ 
perature point between the flow regions for the different 
ensembles. We would like to compare the renormalized 
Polyakov loop obtained with gradient flow to the conven¬ 
tional renormalization of the Polyakov loop based on the 
static potential at zero temperature. We use the contin¬ 
uum extrapolated results for the renormalized Polyakov 
loop obtained using the normalization condition of the 
potential riV{r = ri) = 0.2605 [19]. So we need to match 
the gradient flow scheme with the potential based conven¬ 
tional renormalization scheme. This is done by matching 
the values of the free energy at a single temperature point 
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for Nr = 12 (T « 200 MeV). The other ensembles will 
be shifted by the same amount, which guarantees that 
the cutoff effects from the different Nr are not obscured. 

After performing this shift we show in Fig. the free 
energy of the static charge F 3 = —TlnP^ in the funda¬ 
mental representation for different Nr- At low temper¬ 
atures, T < 200 MeV, we see some Nr dependence and 
the free energy obtained from the gradient flow is larger, 
but approaches the continuum result with increasing Nr- 
The largest deviations from the continuum results are 
about 10% for Nr = 6 and are few percent for Nr = 12. 
For T > 200 MeV the cutoff effects are much smaller 
and we see agreement with free energy obtained using 
conventional renormalization and the results for the free 
energy obtained using gradient flow. Clearly any differ¬ 
ence between the two approaches should vanish in the 
continuum limit. We performed the continuum extrap¬ 
olation of the free energy for different values of the flow 
time and verified that this is indeed the case. In partic¬ 
ular, the deviations between the results obtained in the 
two renormalization schemes that can be seen in Fig. 
disappear after taking the continuum limit. The details 
of this analysis are presented in Appendix where we 
also show explicitly that different choices of the flow time 
amount to a constant shift in F 3 . 

We close this section by noting that we also calculated 
the Polyakov loop in fundamental representation for the 
smaller light quark masses, namely mi = ms/40. Com¬ 
pared to the mi = TOs /20 results we see a downward shift 
of F^. This shift in the free energy is consistent with the 
shift in the deconfining temperature of about 3 MeV, i.e. 
by shifting the TOs/ 40 data by 3 MeV to larger temper¬ 
atures we make them agree with the TOs/ 20 data. The 
mass dependence of F3 is discussed in detail in Appendix 
0 Next we want to use the gradient flow approach to 
calculate the renormalized Polyakov loop in higher rep¬ 
resentations. 


IV. POLYAKOV LOOP IN HIGHER 
REPRESENTATIONS 


We calculated the expectation value of the Polyakov 
loop in higher representations, namely sextet, octet, de- 
cuplet, 15, 15', 24 and 27 using Eqs. ([^-(12) and the 
gradient flow in the same manner as described in the 
previous section. In particular, we used the flow times 
defined by Eq. 0 also here. The numerical results for 
the Polyakov loops in higher representations for Nr = 6 
are shown in terms of the corresponding free energies 
Fn = —T\nPn in Fig. The vertical scale in the figure 
has been shifted by 100 MeV so that the value of F 3 is the 
same as in the previous section. For representation 15' 
the data for flow time 2/o are very noisy around T = 300 
MeV and are therefore not shown in the figure. The free 
energy of the static charge is larger for the higher repre¬ 
sentations at low temperatures. This can be understood 
as follows: At very low temperatures the free energy of 
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Figure 3. The free energy of the static charge in different 
representations as function of the temperature for Nr = 6. 


the static charge is determined by the binding energy 
of the lightest static-light hadron that can screen that 
charge. For the free energy in the fundamental represen¬ 
tation it is given by the mass of the lightest static-light 
meson, for the free energy in the sextet representation it 
is given by the binding of baryon with two static and one 
light quarks, for the free energy in the octet represen¬ 
tation it is determined by the gluelump mass, while for 
higher representation it is determined by binding energies 
of more exotic states. The free energies follow the hierar¬ 
chy that one expects for the hierarchy of the binding en¬ 
ergies of the corresponding hadrons, e.g. the binding en¬ 
ergy of static-light meson can be estimated to be around 
600 MeV, the binding energy of static-light baryon to be 
around 1 GeV [23], while the gluelump binding energy 
is about 2 GeV |25|. The larger values of the free ener¬ 
gies at low temperature lead to smaller signals and thus 
more noisy data. In fact, without the gradient flow it 
is impossible to extract signals for the Polyakov loops in 
higher representations. The gradient flow increases the 
signal by removing the ultraviolet noise if the flow time 
is sufficiently large. Large flow times result in better sig¬ 
nal. Unlike for the fundamental Polyakov loops extract¬ 
ing signals at low temperatures with flow time smaller 
than defined by Eq. 0 is challenging. 

The temperature dependence of the free energy of 
static charges is larger for higher representations as can 
be seen in Fig. At highest temperatures the free en¬ 
ergies in the higher representations are negative and are 
larger in absolute value than the free energy in the fun¬ 
damental representation. This also means that Polyakov 
loops in higher representations are significantly larger 
than one. These features can be understood in term of 
the weak coupling calculations. In leading order pertur¬ 
bation theory the free energy of static charges is 

Fn{T) = -Cnasmo, (18) 

where Cn is the quadratic Casimir operator of represen- 







5 


tation n, as is the coupling constant, and mo ^ y/A^KoTsT 
is the leading order Debye mass. The values of the 
quadratic Casimir operators are given in Table |Tj Accord¬ 
ing to the above equation the free energy of static charges 
satisfies Casimir scaling, i.e. the free energies in vari¬ 
ous representations only differ by the value of C„. This 
Casimir scaling holds in perturbation theory up to or¬ 
der ttg |26l[^ In terms of the Polyakov loops the Casimir 
scaling implies 

Pg = = ... , (19) 

where = Cn/C^. The values of are also given in 
Table H 

Non-perturbatively Casimir scaling of the Polyakov 
loop was studied on the lattice in SU(N) gauge theories 
as well as in two-flavor QCD with heavy quarks. In these 
studies the renormalized Polyakov loop in higher repre¬ 
sentations was calculated assuming the Casimir scaling 
for the renormalization constants of the Polyakov loop, 
i.e. by rescaling the renormalization constants of the fun¬ 
damental representation . This assumption is 

closely related to the Casimir scaling of the zero tem¬ 
perature potentials since the renormalization constants 
are related to the potentials. In SU(3) gauge theory 
the zero temperature potentials between static charges 
in various representations have been calculated |27j . It 
has also been shown that Casimir scaling holds for the 
potentials after subtracting the UV divergent part from 
the potentials to an accuracy better that 5% for dis¬ 
tances r < 1 fm 122!. Furthermore, Casimir scaling of 
the zero temperature potentials holds up to order in 
perturbation theory and its breaking is numerically small 
[28] . The Casimir scaling of the potential in SU(3) gauge 
theory is of course only approximate. For large enough 
distances it is clearly violated since the potential in the 
adjoint and higher representations will saturate at some 
finite value of r due to string breaking, while the potential 


Table I. Values of the quadratic Casimir C„ in representation 
n and its ratio to the value of the fundamental quadratic 
Casimir i?„ = CnjCz- See, e.g., |5] for details. 


n 


Rn 

3 

4/3 

1 

6 

10/3 

5/2 

8 

3 


10 

6 

9/2 

15 

16/3 

4 

15' 

28/3 

7 

24 

25/3 

25/4 

27 

8 

6 


^ To see this it is important to re-exponentiate the perturbative 
expansion of in terms of Fn- 


in the fundamental representation is linearly rising with 
r. In general, Casimir scaling does not hold for the phe¬ 
nomenon of string breaking; string breaking in various 
representations will happen at different distances deter¬ 
mined by the masses of various static hadrons (see the 
discussion above). Strictly speaking the renormalization 
of the Polyakov loop in higher representations would re¬ 
quire calculating and fixing the normalization en(go) of 
the potentials in different representations at zero tem¬ 
perature, i.e. one needs to define separate renormal¬ 
ization constants for each representation independently. 
The choice is just one economical scheme for 

defining the renormalization constants in higher repre¬ 
sentations. 

Using the gradient flow we can calculate the renormal¬ 
ized Polyakov loop in higher representations without any 
assumptions. Furthermore, as already discussed the gra¬ 
dient flow is instrumental for obtaining signals for the 
Polyakov loops in higher representations at low temper¬ 
atures. In fact, we are not aware of any other methods 
that can achieve this in full QCD. 

In Fig. j^we show the renormalized Polyakov loops in 
various representations for AC = 8, 10 and 12 scaled by 
the ratio of the corresponding Casimirs (cf. Eq. @)- 
Here we also impose the additional normalization that 
connects the free energy in the conventional renormal¬ 
ization scheme and the gradient flow renormalization 
scheme. At high temperatures we observe Casimir scal¬ 
ing. This is expected based on the previous lattice stud¬ 
ies. At low temperatures, on the other hand, we see de¬ 
viations from the Casimir scaling. From the above figure 
one can also see that the dependence of the Polyakov 
loop is small. In fact, no Nr dependence is seen within 
errors in Fig. [^ Thus cutoff effects in the Polyakov loop 
in higher representations are under control. 

The deviations from Casimir scaling at low tempera¬ 
tures cannot be seen well in Fig. [^ since the Polyakov 
loops are small there. To clearly see deviations from the 
Casimir scaling at low temperatures we introduce the di¬ 
mensionless combination 

pl/flra 

Sn = l-^. ( 20 ) 

Our numerical results for Sn for various representations 
and various Nr are shown in Fig. [^as a function of the 
temperature. We see that breaking of Casimir scaling for 
T > 250 MeV is of the order of a few percent, but be¬ 
comes significant for lower temperatures. This is the first 
time that breaking of Casimir scaling for the Polyakov 
loop is seen in lattice calculations. In previous studies 
no conclusive statements could be made due to large sta¬ 
tistical errors or large volume effects. The above results 
imply that for T > 250 MeV color screening follows the 
perturbative pattern, while at lower temperatures it is 
strongly nonperturbative. Figure]^ also shows that the 
breaking of Casimir scaling is independent of the value 
of Nr- We demonstrate this in the case of the octet rep¬ 
resentation. This is another way to see that cutoff effects 









6 





Figure 4. Polyakov loops in various representations scaled by the ratio of the appropriate Casimir operators (see text) for 
Nt = 8 (left), Nt = 10 (middle) and Nt = 12 (right). The Polyakov loops have been rescaled by exp(—A 3 /T), A 3 = 100 MeV, 
to match to the conventional scheme for the fundamental Polyakov loop. 




Figure 5. The ratios S„ characterizing the breaking of Casimir scaling for Nt = 6 (left) and various representations n, and for 
the octet representation and various Nt (right). 


in the higher representations are small. 


One may ask to what extent the observed Casimir scal¬ 
ing or its breaking depends on the value of the flow time. 
Different flow time corresponds to different renormaliza¬ 
tion schemes, i.e. to different choices of We cal¬ 
culated Sji for flow times / = /o, 2/o and 3/o and we 
do not see signihcant flow time dependence of this quan¬ 
tity. The corresponding numerical results are presented 
in Appendix Therefore, the above statements about 
the breaking of the Casimir scaling at low temperatures 
are independent on the choice of /. 


We also examined the volume dependence of the 
Polyakov loop in higher representations and did not 
find significant volume dependence. Thus, the observed 
breaking of the Casimir scaling is not a hnite volume ef¬ 
fect. The details of this analysis are given in Appendix 

o 


V. SUMMARY AND DISCUSSION 

We discussed the renormalization of the Polyakov loop 
with the gradient flow. We applied the gradient flow with 
the Symanzik gauge action, i.e. the Symanzik flow, and 
calculated the Polyakov loop after the evolution of the 
gauge fields in flow time up to a fixed value / = 
in physical units, which fixes the renormalization scheme 
for the free energy. With this approach it was possible to 
cover a wide temperature range from temperatures as low 
as T = 116 MeV and up to T = 815 MeV. We compared 
our results for the fundamental Polyakov loop P3 with 
results for the renormalized Polyakov loop obtained in 
the conventional scheme based on the static potential, 
and found very good agreement at all temperatures. 

In addition we calculated the renormalized Polyakov 
loop in higher representations. We found that Casimir 
scaling is approximately fulfilled for full QCD for temper¬ 
atures above T = 250 MeV in agreement with previous 
studies, possibly indicating the weakly coupled nature 
of quark gluon plasma at high temperatures. At lower 
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for the Wilson flow that we used for the Symanzik flow 
(cf. Sec. III). From Fig. [^one finds that at low temper¬ 
atures the cutoff effects are small and both results agree. 
At higher temperatures this changes as for temperatures 
above 350 MeV the Wilson flow produces smaller values 
for the free energies and at some point it is below the 
results obtained using the conventional renormalization 
procedure. This is most likely due to larger cutoff effects 
in the case of the Wilson flow. 


Appendix B: FLOW TIME DEPENDENCE AND 
THE CONTINUUM LIMIT OE THE FREE 
ENERGY 


Figure 6. We compare the renormalized fundamental 
Polyakov loop free energy Fs obtained with Wilson and 
Symanzik flow as described in Sec. |III| In addition we show 
the continuum extrapolated, renormalized free energy from 
|19 | (black triangles). While for lower temperatures the cut¬ 
off effects are small, at higher temperatures the Wilson flow 
shows larger deviations from the conventional result than the 
Symanzik flow. 


temperatures, however, we found for the first time large 
deviations from Casimir scaling. 

The renormalization of the Polyakov loop with gradi¬ 
ent flow is very useful for studying its behavior at high 
temperatures, where performing zero temperature calcu¬ 
lations is very costly. We will discuss this in a forthcom¬ 
ing publication m- 
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Appendix A: COMPARISON OF SYMANZIK 
AND WILSON FLOW 


In this appendix we discuss the flow time dependence 
and the continuum limit of the free energy of fundamental 
charge. As discussed in the main text in the continuum 
limit the free energy of a static charge should be inde¬ 
pendent of the flow time up to an additive temperature 
independent constant, i.e. F„(/) — F’„(/') should be tem¬ 
perature independent. Here we show that this is indeed 
the case using the fundamental free energy as an exam¬ 
ple. To perform the continuum extrapolations we split 
the temperature region used in our study into the low 
temperature region, corresponding to T < 280 MeV, and 
the high temperature region, corresponding to T > 280 
MeV. In these intervals we use the values of / that sat¬ 
isfy the condition given by Eq. (16). We perform inter¬ 
polations of the free energy separately in these intervals 
using smoothing splines and the R package [29]. The er¬ 
rors of the interpolations are estimated by the bootstrap 
method and in some cases adjusted such that they are 
comparable to the statistical errors of the lattice data. 
We perform continuum extrapolations at temperatures 
separated by 5 MeV using the form a + b/N^ + c/N^ and 
the results of these extrapolations are shown in Fig. 

As one can see from the figure, after continuum extrap¬ 
olations F’3(/) — is temperature independent as 

expected. In the figure we also compare our results for 
F3 with the continuum extrapolated results obtained in 
the conventional way. The continuum extrapolated re¬ 
sults obtained for different flow times have been shifted 
by a constant to match the free energy in the conven¬ 
tional scheme. After this shift our results agree with the 
results obtained in the conventional scheme, in particu¬ 
lar there is no discrepancy at low temperature previously 
observed at hxed (cf. Fig. [^. 


In this appendix we compare the renormalized funda¬ 
mental Polyakov loop obtained with the Symanzik flow, 
as described in Sec. jllH with the renormalized funda¬ 
mental Polyakov loop obtained with the same procedure 
but using the Wilson gauge action in Eq. (13), i.e. us¬ 
ing the Wilson flow. In Fig. [^ we show obtained us¬ 
ing Symanzik flow and Wilson flow, and compare these 
with the conventionally renormalized free energy. We 
use the same flow times and same matching procedure 


Appendix C: QUARK MASS DEPENDENCE AND 
VOLUME DEPENDENCE OF THE 
RENORMALIZED POLYAKOV LOOPS 


In this appendix we discuss the quark mass and volume 
dependence of the free energies F„. In addition to the cal¬ 
culations of the Polyakov loop for mi/rris = 1/20 we also 
performed calculations for the smaller light quark mass 
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Figure 7. The difference -F 3 (/) — F^if') for different flow time (upper panels) and the comparison of -F’ 3 (/) with the continuum 
results in the conventional scheme (lower panels). The results in the low and high temperature regions are shown separately in 
the left and right panels, respectively. 


mi/ms = 1/40 on 32^ x 6 and 32^ x 8 lattices. Since for 
Nt = 6 we have two different volumes we can make some 
statements about finite volume effects as well. In Fig. 
we show the temperature dependence of the fundamental 
and adjoint free energies as function of the temperature 
for two different quark masses and flow time 3/o. We see 
that the free energies show some quark mass dependence, 
namely they are smaller for the smaller quark mass. The 
relative difference of the free energies calculated for the 
two quark masses is about the same for triplet and octet 
charges and for Nr = 6 and Nr = 8 . This difference 
may be understood in terms of change in the transition 
temperature. Shifting the ms/20 data by 3 MeV to lower 
temperatures almost eliminates this difference. 

For Nr = 8 the spatial volume is the same for both 
quark masses, but for Nr = 6 the spatial volumes are 
different, namely we use 24^ x 6 and 32^ x 6 volumes. 
Since the shift in the free energies is the same for Nr = 8 
and Nr = 6 and can be understood as a quark mass effect 


we conclude that volume effects in the free energies in the 
fundamental and adjoint representations are smaller than 
the estimated errors and thus can be neglected. Similar 
conclusions can be made for the free energies in other rep¬ 
resentations. Therefore the observed breaking of Casimir 
scaling is not affected by finite volume effects. 

Appendix D: FLOW TIME DEPENDENCE OF 
CASIMIR SCALING 


In this appendix we discuss the Casimir scaling in 
terms of at different flow times. In Fig. |^we show 
Sn in various representations for Nr = 6, 8, 10 and 12. 
We use the values of flow time f = fo, 2/o and 3/o. For 
smaller values of flow time the data are too noisy to al¬ 
low conclusive statements. From the figures we see that 
the flow time dependence of Sn is very small for all Nr- 
Therefore, we conclude that the Casimir scaling or its 
breaking is independent of the flow time. 
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Figure 8. The free energy in the fundamental (triplet) and adjoint (octet) representations calculated for two different quark 
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Figure 9. The measure of the Casimir scaling Sn 
/ = /o, 2/o and 3/o (from left to right). 


shown for N-, 


6, 8, 10 and 12 (from top to bottom) and flow times 
















































































































